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1. Introduction 

One-dimensional diffusion in a potential V(x) is described by the Fokker-Planck equation 

■x) + 2±\f(x)4>(x)] = k*<l>(x) t (1.1) 



where 



dx 2 dx 

/w = -iA ni) . (1 . 2) 

The Fokker-Planck equation (11. ip is equivalent to the Schrodinger equation 
d 2 

ij){x) + V s (x)ip(x) = k 2 ^(x), (1.3) 



dx 2 

where ip and Vs are related to <fi and / in (II. ip by 

0( X ) = e~ v ^{x), V s (x) = f{x) + f\x). (1.4) 

We assume that Imfc > 0. Let Gs(x,y;k) denote the Green function for the 
Schrodinger equation, satisfying 

d 2 



-V s (x) + k 2 



G s (x,y;k) = 6(x-y) (1.5) 



-2—f(x) + k 



G F (x,y;k) = 6(x - y). (1.7) 



<9x 2 

together with the boundary conditions Gg(x,y; k) — > as \x — y\ — > oo for Imk > 0. 
For real k, we define Gs(^, y; k) = lim^o Gg(x, y;k + ie). Without loss of generality we 
assume x > y. The Green function for the Fokker-Planck equation is given by 

G F (x, y; k) = e-^- v ^ 2 G s (x, y; k). (1.6) 

This G F (x,y]k) satisfies 

' d 2 2 d 

dx 2 dx ' 

In previous papers [UI2], we introduced a new method for systematically calculating 
the expansion of Gs(k) in powers of k. In these papers, it was assumed that the potential 
V(x) either steadily tends to a finite limit, or steadily diverges to +oo or — oo, as 
x — > ±oo. We wish to extend our method to be applicable to more general potentials, 
including the cases where V(x) oscillates as x — > ±oo. In the present paper we study 
the simplest of such cases, namely, the cases where the potential is purely periodic. We 
assume that the potential V(x) is a real- valued periodic function with period L, 

V(x + L) = V(x), (1.8) 

and that V(x) is piecewise continuously differentiabl^||. 

The Fokker-Planck equation with a periodic potential has various applications such 
as Josephson junctions [1] and superionic conductors [5j|6], to name a few. And, needless 
to say, the Schrodinger equation with a periodic potential plays an essential role in solid 

X We allow V(x) to have jump discontinuities. Although Vs(x) does not make sense when V(x) is 
discontinuous, the Fokker-Planck equation is well defined even for such V. See [3] or footnote 1 of [2]. 
When Vs does not make sense, Gs is defined by (|1.7j) and (|1.6j) . 
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state physics. Thus, the study of periodic systems is significant in itself. Schrodinger 
operators with periodic potentials have been extensively investigated over the years, and 
much is known about their spectral properties [7-17]. 

To carry out the analysis of the Green function, we adopt a formalism that 
deals with reflection coefficients. Reflection coefficients are of extreme importance in 
scattering theory, and their properties have been studied by many researchers, especially 
in connection with the inverse scattering method [12-16]. The formula derived in [TJ 
provides a new method for calculating the expansion of reflection coefficients up to an 
arbitrary order in k. The aim of the present paper is to modify this formula and make 
it applicable to the periodic case. The resulting formula serves as a new tool for the 
analysis of periodic systems, and enables us to deal with periodic and non-periodic 
potentials on the same basis. Although we restrict ourselves to periodic potentials in 
this paper, this method can be further extended to more general potentials which show 
oscillatory behavior at infinity. 

The necessary results from the previous papers will be briefly reviewed in section 3, 
after making some definitions in the next section. In sections 4-8, we discuss how to 
adapt the expansion formula to the periodic case. Using the expansion of reflection 
coefficients, we shall derive the low-energy expansion of the Green function in sections 9 
and 10. As is well known, the energy spectrum of a periodic system has a band structure. 
The expansion in powers of k is an expansion from the bottom of the lowest band, and 
it is effective for calculating the Green function in the lowest band. An example is given 
in section 11. 



2. Preliminaries 



In this section we summarize the necessary definitions and notation. We define the 
evolution matrix U(x, x'; k) as the 2x2 matrix satisfying the differential equation [17] 

±UM *)=("* '£> (*,*■; k) (2.1) 
with the initial condition 



U(x' } x';k)=[ " . (2.2) 




It is obvious from (12. ip that the upper-right and lower-right elements of U are obtained, 
respectively, from the lower-left and upper-left elements by changing the sign of k. We 
write these elements as 

U(x, x'; k)=( a ^ x '' k ) P(x,xf;-k)\ (23) 
I /3(x,x';k) oc{x, x'; — k) J 

From ( 12. Ill . we can easily show that a(x, x'\ k)+(3(x, x'\ k) and a(x, x'\ —k)+/3(x, x'\ —k), 
as functions of x, are solutions of the Schrodinger equation ( II. 3p . The matrix equation 
( 12. ip is thus equivalent to ( 11. lj) and ( 11. 3p . 
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Since the matrix multiplying U on the right-hand side of (12.1 j) is traceless, the 
determinant of U is unity. Namely, 

a(x, x'\ k)a(x, x'; —k) — (3(x, x'\ k)(3(x, x; —k) = 1. (2.4) 
It is easy to see that the evolution matrix has the property 

U (x 3 , x 2 ; k)U(x 2 , x 1 ;k) = U [x 3 , x u ; k). (2.5) 
If the potential is a periodic function satisfying ( 11 .81) . then f(x + L) = f(x), and so 

U{x + L,x' + L;k) = U{x,x';k). (2.6) 

The transmission coefficient r, the right reflection coefficient i? r , and the left 
reflection coefficient Ri are defined in terms of a and /3 as 

r(x, a/; fc) = - 1 , (2.7a) 

R r (x,x;k) = — — , Ri{x,x;k) = — — . (2.7b) 

a(x,x';k) a(x,x';k) 

(See [H]. Similar construction of the reflection coefficients is used, for example, in [IB].) 
Let us explain why these quantities are called the transmission and reflection coefficients. 
We consider a potential which is identical to V(x) inside the interval (xi, x 2 ) and defined 
to be constant outside this interval. Namely, we define 

{V(x\) (x < x\) 
V(x) (x 1 <x<x 2 ) f XuX2 (x) = --—V XuX2 (x). (2.8) 
V(x 2 ) (x 2 < x), 

The Schrodinger equation corresponding to this Fokker-Planck potential V Xl>X2 is 

- ^ + (fl*. + fLJ * = (2-9) 

As shown in appendix A, this Schrodinger equation has two independent solutions that 
behave outside the interval ( 

. r{x 2 ,x 1 ;k)e- ik( - x - x ^ (x < x x ) 

ipi(X) — < ., . , ... . (2.10a) 

e- lk ( x ~ X2 ^ + R r {x 2 ,x 1 ;k)e lk( - X - X2 ) (x>x 2 ), 

e *C*-»0 + R t (x 2 , Xl ; k) e~ ik ^ (x < x x ) 

ril> ' T(x 2 ,Xi,k)^ X - X ^ (X>X 2 ). { ' ' 

Thus, R r (x 2 ,xi, k) (or Ri(x 2 ,xi, k)) is the factor multiplying the wave reflected from 
the interval (x\,x 2 ) when an incident wave is coming into this interval from the right 
(or, respectively, left), and t(x 2 ,xi, k) is the factor multiplying the transmitted wave. 

§ In some of the previous papers, there is an error in the equations corresponding to ()2.10|) (equations 
(2.3) of [2], equations (1.6) of p], and equations (3.2) of 20 ). The functions defined in these 
equations are solutions of the Schrodinger equation, and not the Fokker-Planck equation. For these 
functions to become solutions of the Fokker-Planck equation, there should be a factor e^ v ^ X2 ' v ^ Xl ^^ 2 
or e~^ v ^ x ' 2 )~ v ^ Xl )y 2 in front of the transmission coefficient t(x2,£i; k). This error does not affect any 
of the results of these papers. 
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We define 

Qf 7 \ R r (x,-oc;k) Ri{oo;x;k) 

Sr{X, k) = — —. jT, Si{X, k) = — rr, 2.11 

1 + Rr(x, — oo; k) 1 + Ri(oo, x; k) 

S(x,k) = S r (x,k) + Si(x,k). (2.12) 

The reflection coefficients for semi- infinite intervals, which appear in (I2.1ip . are well 
defined as long as Imk > 0. (The well-definedness of R r (x,—oc;k) for periodic 
potentials is explained in appendix B. For non-periodic potentials, see appendix G 
of j2]. A related discussion is found in [19].) For Irak = 0, we define them as 
R r (x, — oo; k) = lini4o R r (x, — oo; k + ie) and i?z(oo, x; k) = lim^o Ri{oo, x;k + ie). As 
explained in appendix C, the functions S r , Si and S are closely related to the Weyl- 
Titchmarsh m-function. The Green function can be expressed in terms of this S as 

m\ 



ik(x — y) — ik I S(z,k)dz 



(2.13) 



G$(x, y: k) = ; exp 

2iky/[l -S{x,kW-S(y,k)] 

In our formalism, we deal with the scattering coefficients r, R r and Ri in an 
generalized form with an additional variable W. First, we define a and (3 by 

a(x,x';W;k) /3(x,x';W;—k) 
j3(x, x'\ W; k) a(x,x';W;—k) 



cosh — — sinh w \ ( a(x,x';k) /3(x,x';—k) 

- sinh E^M cosh E^pEl )\p( x , x '-k) a(x,x';-k) 



(2.14) 



(Here, as elsewhere in this paper, the bar does not denote complex conjugation.) Then, 
f, R r and Ri are defined in the same way as (12. 7p with a — > at and — > j3: 

( >wir\= 1 _ y/l-e(x,W)r(x,x';k) 

T[X,X ' ' a(x,x';W;k) ~ 1 - £{x, W)R r (x, x>; k) ' 1 ^ 

R (x x'- W k) = - Rr(x,x';k)-gx,W) 

Kr[X,X,W,k) - a(x,x';W;k) ~ l-£(x, W)Rr(x, x'; k) ' [ZAbt>) 

^/ i t t 7" 7 \ B(x,x':W;-k) n , . f (x, W)t 2 (x, x'\ k) . . 

Ri(x, x'; W; k) = - _V \ ' J = Ri(x, x'; k) + - q \ ±-±± 2.15c 
a[x, x'\ W ; k) 1 — ^[x, W )Rr{x, x , k) 



W - V(x) 

£(x, W) = tanh { ] . (2.16) 



where 



The meaning of these generalized scattering coefficients is explained in [2T]. The original 
scattering coefficients r, R r and Ri are recovered from f, R r and Ri by setting W = V(x). 
For k = 0, equation (12.11) can be exactly solved. It is easy to see that 

a(x,x';0) = cosh V ^~ V ( X \ 3{x, x'; 0) = sinh Yi*lzX&l . ( 2 .17) 
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Hence we have the expressions of f, R r and Ri for k = as 

W - V(x') 

f(x, x'\ W; 0) = sech 2^ 

W - V(x') - W - V(x') 
R r (x, x'; W; 0) = - tanh K -^, x'; W; 0) = tanh K -^. (2.18) 



For n = 1, 2, 3, . . . and —00 < a < b < 00, we define the notation 



[cxi, cr 2 , . . . , a n ] b a = ■■ dz 1 dz 2 ■ ■ ■ dz n exp ^ crjV^ 

J J a<z\<Z2<--<z n <b *-j=l 

where each <7j is either +1 or — 1. For simplicity, we use the obvious abbreviation [+} b a , 

Ha, [H — la) etc in P lace of [ +1 ]a> [~l]o» _1 la> etc - Tne integrals of the form of 
( I2.19P satisfy the multiplication rule 

[(Ti, cr 2 , . . . , a n } b a x [er'] = [er', cr 1; cr 2 , . . . , cr n ] a + [ai, a', cr 2 , . . . , a n } b a + [cxi, o 2 , er' . . . , a n } b a 

+ ■■■ + [at, a 2 ,..., a', a n ] b a + [a x , a 2 , . . . , a n , a'} b a . (2.20) 

When the potential satisfies (II. 81) . the quantities defined by 

f e~ v ^dz, P ee [+f x _ L = f e v ^dz (2.21) 

Jx—L J x~L 

are independent of x. It is also convenient to define 

L = VPM, e Vo = y/P/M, i.e. V Q = (1/2) log(P/M). (2.22) 

From f l2~20|) and f l2~2T|) . we have 

= ^/2, = M 2 /2, + [-+) x x -l = PM = Lq. (2.23) 

We can check that the left-hand sides of ( I2.23P are independent of x by differentiating 
them with respect to x, using 

A [ au a2 , . . . , = [a u a 2 ,..., a n _^_ L e*» v <*> - [<t 2 , a,,..., a n \l_ L e^H (2.24) 

We define the operators A and B, which act on functions of x and W, as 

Ah(x, W) = -0-h(x, W), (2.25) 
ox 

Bh{x, W) = ^- {sinh[H/ - V(x)]h(x, W)} 

cosh[H/ - V(x)} + smh[W - V{x)\— \ h(x, W). (2.26) 

We assume that the functions h(x, W), on which these operators act, are (i) piecewise 
continuously differentiabl^jj] with respect to x, and (ii) analytic with respect to W on 
the real axis. 

Finally, the operator T> is defined as 

Vh(x,W) = h(x,W) -h(x,V ), (2.27) 

Here we do not require h(x, W) to be continuous in x. We allow h to have jumps, and so dh/dx 
may contain delta functions. The argument in pQ is valid without the requirement of continuity. 



^3) 



(2.19) 



M = [-]*_ L 
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with Vq defined by (I2.22p . In terms of differential and integral operators, we can write 

rW q 

V = k iW W' <2 ' 28) 

3. The low-energy expansion formula for non-periodic potentials 

Here we review the expansion formula for reflection coefficients derived in pQ , where the 
potential V(x) was assumed to be monotone for sufficiently large \x\. For details, see [1] 
and [2]. In this section, V(x) is not assumed to be periodic. 

With expression (12.131) . the analysis of the Green function is reduced to the 
analysis of the reflection coefficients R r (x, — oo; k) and Ri(oo,x;k). Here we consider 
R r (x,—oc;k). As a function of x, the reflection coefficient satisfies a nonlinear 
differential equation of Riccati type. By introducing the generalized reflection coefficient 
( 12.1551) with an additional variable W, we can turn this nonlinear equation into a linear 
partial differential equation for the two variables x and W. This equation has the form 

(A - 2ikB) [Rr{x, -oo; W: k) + W)] = [1 - ^(x, W)}f(x) (3.1) 

with the operators A and B defined by (12.251) and (12.261) . respectively. The reflection 
coefficient can be obtained by solving (13. ip . Solving the differential equation (13.11) with 
an appropriate boundary condition is equivalent to finding the inverse of the operator 
A — ikB with an appropriately restricted domain. 

Let us start with the case k = 0. If we restrict the domain of the operator A to 
functions satisfying 

lim h(x,W) = 0, (3.2) 

X — OO 

then we have the inverse of A as 

A- 1 g(x,W)= [ g(z,W)dz, (3.3) 



where g(x, W) = (d/dx)h(x, W) with some h satisfying (I3.2p . The operator A^ 1 given 
by (13.31) satisfies A~ 1 Ah = h and AA~ x g = g if h satisfies (13. 2p . 

This can be extended to k ^ 0. We restrict the domain of the operator A — ikB to 
functions h(x, W) which satisfy 



-oo 



for any x, where 



-/ u/ n - t/^ , l ^ + Ri(x,z;W;k) 

u(x, zs W- k) = V{z) + log Y^^^y (3-5) 



Then the inverse of A — ikB is given by 

-X -2 



dz l L; J.M z, Q(x, z- W; k)). (3.6) 

.„„ I - Rf(x,z;W;k) 

This (A - ikB)- 1 satisfies (A - ikB)-\A - ikB)h = h and (A - ikB) (A - ikB)- x g = g 
provided that h satisfies (13 .4p . In (I3.4p . the expression h(z, cu(x, z; W; k)) stands for the 
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quantity obtained from h(x, W) by inserting z and Cu(x, z; W; k) in place of x and W, 
respectively. From (12 . 1 8[) and (13.51) . we have u(x, z\ W; 0) = W. By using this and 
( 12.18!) . we can see that (13-41) reduces to (13.21) if we set k = 0. Similarly, the right-hand 
side of (I3.6P reduces to (13. 3p when k is set to be zero. 

It can be shown that R r (x, — oo; W; k)+^(x, W), as a function of x and W, satisfies 
(13. 4p . So it is possible to apply (^4 — 2ikB)~ l to both sides of (13.11) and obtain 

R r (x, -oo; W; k) = -£ + (A- 2ikB)- 1 (l - £ 2 )/. (3.7) 
By iterating the identity 

(A - 2ikB)- 1 = A' 1 + 2ik{A - 2ikB)- l BA-\ (3.8) 
we have the expression 

(A - 2ikB)~ 1 = [1 + ikC + (ik) 2 C 2 + ■■■ + (ik) N C N ] A' 1 

+ (ik) N+1 (A - 2ikB)- 1 AC N+1 A-\ (3.9) 

where 

C = 2A~ 1 B. (3.10) 

The expansion of R r (x, — oo; W; k) in powers of k can be obtained by substituting ( 13.91 ) 
into (15771) . We have 

R r = f„ + i/cn + (iA;) 2 f 2 + • • • + (i/c^iv + Ptv (3.11) 

with the coefficients of the expansion 

f = ^- 1 (l-e 2 )/-^, (3.12a) 

f n = C n (f + (n>l), (3.126) 

and the remainder term 

= (ik) N+1 {A - 2ikB)- 1 Ar N+1 . (3.13) 

This method works well if V(x) either converges to a definite value or diverges to 
±oo as x — > — oo. From ( 11.21) and ( 12.161) we find 

Substituting this into ( 13. 3p gives 

^d-a/=r fi tot ff - y wu ={ -^ ff - y '-°°'. (3.i5) 



dz 2 J 2 

The zeroth-order coefficient of the expansion (equation (13 . 1 2 ) is thus 

f = -tanh ^ '-. (3.16) 

If V(— oo) = ±oo, then f = ±1. The higher order coefficients can be obtained by 
substituting (I3.16P into (j3.12fr|) . The resulting expressions are given in [TJ and [2]. 
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4. Application of the expansion formula to periodic potentials 

The expansion formula shown in the previous section cannot be directly applied to 
periodic potentials. If V is periodic, the right-hand side of (13. 16[) does not have a 
definite value. It is the objective of this paper to make a modification to this formula 
so that it becomes applicable to the periodic case. 

Equation (13. ip holds for any potential, whether it is periodic or not. The function 
h(x, W) = R r (x, — oo; W; k) + W) is periodic in x when the potential is periodic. If 
h is periodic, then (13.41) is satisfied as long as Imk > 0. Therefore, the basic expressions 
( 13. TP and (13.61) are correct even when the potential is periodic. (When Imk = 0, it is 
necessary to replace k in the right-hand side of (13. 6p by k + ie with e > 0, and take the 
limit e — > after evaluating the integral.) 

If we try to apply the method of section 3 to a periodic potential, we need to let 
A -1 act on periodic functions, as in equations (I3.12al) and (13. 12 £>[) . However, since (13.21) 
is not satisfied when h(x, W) is periodic, we cannot use (13.31) for periodic functions. The 
right-hand side of (13. 3p is indeterminate when g(x, W) is periodic in x. 

Thus, in order to deal with periodic potentials, we must discard (13 .31) and find 
another appropriate expression for the inverse of A. Since (13.61) is valid for both periodic 
and non-periodic cases, we can derive A" 1 from (13.61) as 

A~ x g{x, W) = \im(A - 2ikB)~ 1 g(x, W) 

fc->0 

As we shall see later, this gives the correct inverse of A even when the potential is 
periodic. This expression includes ( 13. 3 p as a special case. As noted below equation 
( 13. 6p . the right-hand side of ( 13. 3 p is obtained by letting k — > inside the integral, i.e. 

" x t 2 (t 7- W- k) r x 

lim ^ ' ' '7 - g(z, u{x, z; W; k)) dz = / g{z,W)dz. (4.2) 



k-*ol-R*(x } z]W;k)' 

Hence, we find that ( 14. ip reduces to (13.31) if the limit and the integral in (14. lj) are 
interchangeable. This is what happens in the non-periodic cases discussed in the 
previous papers. In the periodic case, however, the limit and the integral cannot be 
interchanged, and so the correct inverse of A is different from ( 13.31) . 

What we mean by "correct inverse of A" is the operator A" 1 that can be used in 
equations (13. 121) to produce the correct expansion of R r . As we have not yet specified 
the domain of A, the term "inverse of A" does not yet have a definite meaning. For the 
time being, let us accept ( 14.11) as the definition of A^ 1 . It will be shown in section 6, 
after properly defining the domain of A, that this A' 1 is indeed the inverse of A. 

The method of section 3 remains valid for periodic potentials if ( 13. 3p is is replaced 
by ( 14. ip . Our main task is to calculate the right-hand side of ( 14.11) to derive an expression 
of A^ 1 for the periodic case. This will be done in the next section. The result is 

A- x g(x,W)= - . * ... V [[ dzdz'smh[V(z')-W]g(z,W), (4.3) 

n smh.[VV — v ) J Jx-l<z<z'<x 
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where L$, Vq and T> are defined by ( I2.22p and (I2.27p . More explicitly, this means 

A' x g{x, W) = - , / dz / dz'\ smh[V(z>) - W]g(z, W) 

L smh(W - Vq) J x _ l J z I 

-smh[V(z')-Vo}g(z,V )y (4.4) 

This, instead of (13.31) . is the inverse of A appropriate for periodic potentials. 

Hereafter we shall always assume that V(x) is periodic. The symbol A" 1 shall 
always refer to (14. ip or (14. 3p . and not (13.31) . 

5. Derivation of (14.3ft 

In this section, we shall derive (14. 31) from (I4.ip . We assume that g(x, W) = Ah(x, W) 
with some function h satisfying h(x + L, W) = h(x, W). Then it follows that 

g(x + L,W) = g(x,W), (5.1) 



g(z,W)dz = 0, (5.2) 

I x—L 

where x is arbitrary. As it turns out, ( 14. ip does not make sense unless (15.2p is satisfied. 
The expressions appearing in ( 13 .6p and (13. 5p can be written as 

f 2 (x, z; W; k) 1 



1 - Rf(x, z; W; k) a 2 (x, z; W: k) - f3 2 (x, z; W; -k) ' 
1 + Ri(x, z\ W; k) a(x, z; W; k) - (3(x, z; W; -k) 



(5.3a) 
(5.35) 



1 - Ri(x, z; W; k) a(x, z; W; k) + (3(x, z; W; -k) 

So, in order to calculate the right-hand side of (14. ip . we need to know the small- 
k behavior of a(x,z;W;k) and j3(x, z; W; — k). By definition, a(x, z;W; k) and 
f3(x, z; W; —k) are obtained from the elements of the matrix U(x, z; k). 

We divide the domain of the integral in ( 14. ip into regions of length L as 

°° px—nL 

= E / • < 5 - 4 > 

> n =0 J x-(n+l)L 

To deal with each term on the right-hand side of ( 15 .4p . we assume that 

x — (n + 1)L < z < x — nL. (5.5) 
Using ( 12. 5p and (12. 6p . the matrix U(x, z\ k) is factorized as 

U(x, z; k) = U(x, x — nL; k)U(x — nL, z; k) = [U(x, x — L; k)] n U (x — nL, z; k). (5.6) 

We can calculate [U(x, x—L; k)] n , which appears in (I5.6p . by diagonalizing U(x, x — L; k). 
The eigenvalues of U(x,x — L; k) are A(A;) and 1/X(k), where 

X(k)=Y(k)-Wl-[Y(k)} 2 , (5.7) 
Y{k) = -[a(x, x — L; k) + a{x, x — L; —k)\. (5.8) 



[U(x,x-L)] r 
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We used ( 12. 4p to derive (15. 7p . As shown in appendix D, the quantity Y(k) is independent 
of x, and so is X(k). The branch of the square root in (15.71) is chosen so that \X(k)\ > 1 
for ImA; > 0. The two eigenvalues of U are the reciprocal of each other, reflecting 
the fact that the determinant of U is unity. Also note that A(— k) = 1/X(k). It is 
straightforward to show that [U(x,x — L; k)] n is expressed in terms of A as 

A~ n (A - a) - X n {X- 1 - a) (3'{X n - X~ n ) 

A - A- 1 I (5{X n - X~ n ) X n (X -a)- A^A" 1 - a) 

(5.9) 

where we have omitted the dependence on k, and used the abbreviation 

a = a(x, x — L; k), (3 = (3(x,x — L; k), /3' = /3(x, x — L; — k). (5.10) 

In general, for any finite X\ and x 2 , we can expand a(x 2 ,xi, k) and f3(x 2 ,Xi] k) in 
powers of k as 

a{x 2 , x 1 ;k) = cosh ^™ - | ( e ^ +v ^ 2 [-]£ + e -^ +v ^ 2 [+]**) 

JL ( e WW-VtoW + e lV(z 2 )-V( Xl )]/2 + 0( p ); (511a) 

(3(x 2 , x i; k) = sinh v ^)- v &) + | ( e ^(*0+v(* a )]/2 [_]*, _ e -[V(*i)+v (sea )]/2 [+ Q 

( e [V(xO-V(* 3 )]/2 [_ + ]« _ e [v(x 3 )-v (a!1 )]/2 [+_]**) + o(P). (5.116) 

The derivation of (15.111) is given in appendix B. With x 2 = x and x\ = x — L, equation 
(15.1 lap reads 

a{x, x - L; k) = 1 - - (e v{x) M + e~ v{x) P) - —PM + 0{k 3 ) 

i.2 

= 1 - ikL cosh[y(x) - V } - —L\ + 0(k 3 ), (5.12) 
where we have used (12T221 and (l2"^3"p . Substituting (15U21 into (EHJ) and (jOl . we find 

K(A;) = 1 - yL 2 + 0(A; 4 ), V / HW = kL + 0(k 3 ), (5.13) 

X(k) = l-ikL ~^-Ll + O(k 3 ), [X(k)}- 1 = l + YkLo-^-Ll + 0{k 3 ). (5.14) 

Hence, the small-fc expressions of the quantities appearing in (15. 9p are 

A - A -1 = -2ikL + 0(k 3 ), (5.15a) 

A — a(x,x — L; k) = 2ifcL sinh 2 V ° + o(k 3 ), (5.156) 

A" 1 - a(x, x-L-k) = 2ikL cosh 2 V ° + o(& 3 ). ( 5 .i 5c ) 

The expressions for /3 and are given by ( 15. 11 6ft as 

0{x, x-L- ±k) = ±ikL smh[V(x) - V ] - y ( [-+]*_ £ - + 0(A; 3 ). (5.16) 
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By substituting (I5.15P and (I5.16P into (15. 9p . we obtain the elements of U(x,x — nL; k) = 
[U(x, x — L; k)] n as 

a(x, x - nL; ±k) = cosh 2 A ±n - sinh 2 X Tn + 0(k 2 ), (5.17a) 

(3(x, x - nL; ±k) = =p cosh ^h^ 1 sinh ( A ™ " A ~") 

- ( [-+H-L - [+-} X X - L ) (a- - a-") + o(e). (5.m) 

(In these equations, A stands for X(k) and not X(±k). Recall that X(—k) = 1/X(k).) 
Since we are going to take the sum over n as in (15 .4p . we cannot expand A™ and A _n in 
powers of k and neglect the higher order terms. Even if \k\ is small, the higher order 
terms are not negligible when n is large. So, we leave X n and X~ n as they are. 
Up to order k, the elements of U(x — nL, z; k) are given by (15.1 1 j) as 

a(x - nL, z; ±k) = cosh YMtpH T | ^v(x )+ v( z )]/2 + e -\y&+v(*)]/2 

+ 0(k 2 ), (5.18a) 

f3(x - nL, z; ±k) = sinh ± [ A ( e W*)+vW]/2 _ e -[v( x )+v(z)]/2 

+ 0(A; 2 ). (5.186) 

The expressions of a(x, z; ±fc) and (3(x, z; ±fc) are obtained by substituting (15. 17ft and 
(I5.18P into U (x, z; k) = U(x, x — nL; k)U(x — nL, z; k), or, explicitly, 

a(x, z; ±k) = a(x, x — nL; ±k)a(x — nL, z; ±k) + j3(x, x — nL; ^fk)f3(x — nL, z; ±fc), 
f3(x, z; ±k) = j3(x, x — nL; ±k)a(x — nL, z; ±k) + a(x, x — nL; ^fk)(3(x — nL, z; ±fc). 

(5.19) 

Substituting the resulting expressions into the definition of a and j3 (equations ( 12.14(1 ) 
yields 

a(x, z; W; k) = (cosh cosh ^) A ™ ~ ( sinh sinh ^) A ~" 

JA [ cosh K^Vo ( e [V( Z ) + V ]/2 { _ ]: -nL + e -[VW+Vo]/2[ + ]*-^ ^ 



+ sinh W ~ V ° ( e lV(z)+V ]/2^x-nL _ e -\V(z)+V ]/2^x-nL\ ^-n 



ik 



a 



sinh ( [-+r x _ L - [+-)l_ L ) (X n - X- n ) + 0(k 2 ), (5.20 

4_Lo 

(3(x, z; W; -k) = (sinh cosh ^r*) A ™ ~ ( cosh sinh 

_^ Ugh ^[V^+Vo^x-nL _ e -[KW+Vo]/2 [+ ]*-^ ^ 

4- sinh ^-^o ( e [V(*)+K.]/2[_]x-nL _j_ g-|V(*)+Vo]/2[ + ]s-n2A ^-n" 

,/,; cosh ( [-+]*_ L - [+-]»_ J (A" - A-«) + 0(A; 2 ). (5.206) 



4L 
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Details of the calculation are given in appendix E. From (15. 3p and (I5.20p we obtain 



r 



f 2 (x,z;W;k) 1 

1 - Rf(x, z; W; k) ~ cosh 2 1 - ( Co7 ™)2 



x j | +ife l±£22! e - v °[+lr" L + ife T . C0l \ v °[-]T nL + kB x 



col' 



1 + c 7 r 



l + ift(a;,*; W;fc) 
1 — Ri(x, z; W; k) 



+ 0(k 2 



-V(z)+V _ 



(5.21a) 



c 7 



1 + CoY' 



+ 0(A; 2 ), 



1 r 

V 0][ + ] x-nL + iA .__^J_ e Vo[_]x-nL + ^ 



1 + c 7 r 



where B\ and I?2 are quantities independent of z, and 



c 



tanh • 



1=^2 =l + 2ikL + O(k 2 ). 



(5.216) 



(5.22) 



As mentioned before, we cannot let 7" ~ 1 + 2nikLo in (I5.2ip . We must take the sum 
over n before using 7 ~ 1 + 2ikL . The z-independent terms kBi and kB 2 in (I5.2ip 
derive from the terms involving ([ — \-}^_ L — [H — ]x-l) * n Q5.20j) . As we shall see, they 
are irrelevant to the final result. 

We assume that g(x, W) can be expanded in powers of e w as 



g(x,W) = Y,9i(x> jW - 

3=0 

Condition (I5.2p requires that 



gj(z) dz = 0. 



(5.23) 



(5.24) 



x—L 



For each term in (I5.23p . expression ( 13. 6 j) gives 
(^-2iA;i3)- 1 ^(x)e jW 



r 2 (x,z;iy;fc) , y(g) 
1 -i? 2 (x,z;iy ; A;) 



l- Rl (x,z;W;k)\ 9 > {Z) - (5 ' 25) 
We divide the integral as (15. 4p . and shift the variable of integration by z — > z — nL. 
Thus, the right-hand side of (15. 25ft is rewritten as 

00 px—nL 00 px 

C(z)dz = J2 C(z)dz = J2 C(z-nL)dz, (5.26) 

n=0 Jx-{n+l)L n=0 Jx-L 

where C(z) stands for the integrand in (I5.25p . For each n in the sum of ( I5.26p . we can 
use (I5^TD . Substituting fl5^il . and using = [+]*, {-f'-nL = Hz, we obtain 

(A- 2ikB)- 1 g j (x)e iW 



JV 



cosh 2 M 

2. n=0 



7" 



X - L 1 - (c 7 



n\2 



1 " c 7 r 
1 + c 7 r 



n \ 



l + ik(l-j)l±^e- v °\+r 



1 r ~,n 

+ ik(l + j)- °—e Vo [-f z + kB 3 



1 - c 7 
2 



#(z) dz + 0(£; 



(5.27) 
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where B 3 (which derives from B\ and B 2 of (I5.2ip ) is independent of z. On the right- 
hand side, the term of order k° vanishes on account of (I5.24p . The part involving B 3 
vanishes for the same reason. Therefore, (15.27!) becomes 

(A- 2ikB)- 1 g j (x)e jW 



ik 



cosh 2 ^ n 

2 n=0 



3i 



r 



(co7 



n\2 



- v - i) 



1 



U + 1) 



1 



*7> 



1 + c 7 n 
i-i 

x-L 



3+1 



Hz9j( Z ) dz 



x—L 



[+) x z 9j(z)dz 



+ 0{k)\. (5.2 



. 1 + c 7 n 

We have |7| < 1 as long as Imfc > 0. (Recall that R r (k) = lim^o-Rr^ + ie) for 
real k.) For an arbitrary integer m > 0, we can calculate the infinite sum 

ik 1 



n=0 



n(l+m) 



7 



l+m 



2(m + l)L 



+ 0(k), 



where we have used 7 = 1 + 2iA;L + 0{k 2 ). Hence, for any analytic function h, 

-1 



lim ik l n hh r 



n=0 



2L 



h(x) dx, 



Jo 



(5.29) 



(5.30) 



as can be proved by Taylor expanding h(x) and using (I5.29P for each term. We can take 
the limit k — > of f !5.28|) by using (I5.30p . To calculate the integral, we use the formula 



1 — CqX 



1 - (c x) 2 \1 + C X 



dx 



2mcn 





\ mi 
- c \ 









e m(W-V ) _ l 

2mtanh^^ 



(5.31) 



From fl538|) . flCTj) and f l53Tj) we obtain 



lim(A-2ikB)- 1 g 3 (x)e ]W 



2L sinh(H/-y ) 



(eO-Dw _ e O--DVo) / [+]^.( z )dz 



x—L 



1 



2L sinh(H/-y ; j.-z 
Taking the sum over j, we arrive at the resul 

1 



(e~ w i+r z 



[-)z9j(z) dz 

x—L 

) 9j (z)^ w dz. 



A~ x g{x,W) 



-V 



x—L 



(e- w [+r z 



-]*)g{z,W) dz. 



(5.32) 



(5.33) 



2L sinh(W -V ) 

Substituting [+]* = f* e v W dz' and [-]* = f* e~ v(z,) dz' into (jOSjl gives (Ojl . 

The above derivation shows that linifc^o(-4 — 2ikB)~ 1 g makes sense and is equal to 
the right-hand side of (14. 3 p if both ( 15. ip and (15. 2p are satisfied. If g satisfies only (15. ip 

If Here we are assuming that A~ x J2j 9j^ W = A~ 1 gje' w . If we can write g(x, W) — Ah(x, W) 
with h(x, W) = J2j hj{x)eP w , this assumption is equivalent to (d/dx) J7 h^ w = J2j(9hj/dx)e : ' w . It 
is not necessary to verify this assumption here, as it will be checked later (in section 6 and appendix F) 
that (|5.33p is correct as long as g(x, W) belongs to the range of A defined in the next section. 
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and not (15. 2p . then (A — 2ikB)~ 1 g contains a term of order k~ x (see (I5.27P and (joTc 
and so it is impossible to take the limit k — > 0. 

Since (d/dx)A~ 1 g = g (see section 6 and appendix F), we can express A~ l g as 

A~ 1 g(x,W)= I g(z,W)dz + c(x ,W), (5.34) 

JXQ 

where x is an arbitrary number, and c(x , W) is independent of x. Setting x = x in 
(I5.34p . we find c(xo, W) = (A~ 1 g)(xo, W). Equation (15.341) can also be directly derived 
by the method described in this section if we start with f x = f x + V [ x ° J 1 ,-,, T instead 

J J —oo Jxq ' Jxq — (n+l)L 

of dS3D. 

6. Domains of A and ^4 1 

Let h and g be functions such that g(x, W) = Ah(x, W). We are assuming that h(x, W) 
is piecewise continuously differentiable in x and analytic in W on the real axis, and also 
that 

h(x + L,W) = h(x,W). (6.1) 

As shown in the last section, A _1 g(x, W) makes sense if (15.11) and (15. 2p are satisfied. 
Equations (15.1 ft and (15.21) follow from ( 16. ip . But this is not enough. For (14.31) to give 
the proper inverse of A, the domain of A must be further restricted. 

As an inverse, A~ x is required to satisfy AA~ x g = g and A~ x Ah = h. From the 
first line of ( 14. ip . it follows that 

lim(^ - 2ikB){A - 2ikB)- 1 = AA- 1 - 2i lim kBA~\ (6.2a) 

k — ^0 k — ^0 

lim(A - 2ikB)- 1 (A - 2\kB) = A' 1 A - 2i lim kA~ l B. (6.26) 

The left-hand sides of (16.21) are identity operators. Therefore, AA~ l g = g holds if 
BA~ l g makes sense, and A~ l Ah = h holds if A~ x Bh makes sense. The former condition 
is always satisfied since A~ x g{x, W) is analytic in W. So, AA~ l g = g holds without 
any further assumption on g. 

On the other hand, the condition for A~ l Ah = h imposes an additional restriction 
on h. We already know that A~ x g makes sense if (I5.2p is satisfied. Replacing g by Bh, 
we can see that A~ l Bh makes sense if 



Bh(z, W) dz = 0. (6.3) 

' x—L 

Namely, A^Ah = h holds if h satisfies ( 16.31) . 

In appendix F, it is verified by direct calculation that the operator A^ 1 given by 
(14. 3 p (or (I5.33P ) satisfies A' 1 Ah = h and AA^g = g, provided that h and g satisfy the 
conditions specified above. Hence, we know that ( 14.31) is indeed the correct expression 
for the inverse of A. 

The conclusion of this section is as follows. We take the domain of the differential 
operator A = d/dx to be the set of functions h(x, W) which satisfy (16.11) and (16.31) . and 
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which are piecewise continuously different iable in x and analytic in W on the real axis. 
Then (14.31) gives the inverse of A, satisfying A~ x Ah = h. The domain of A" 1 is the 
range of A. It consists of functions g(x, W) which satisfy ( 15. ip and ( 15.21) . and which 
are analytic in W and continuous in x except possibly for some finite jumps and delta 
function singularities. 

Let us note that R r (x, — oo; W; k) + £(x, W), which appears on the left-hand side 
of (I3.ip . belongs to the domain of A defined above. It is obvious that R r + £ satisfies 
( 16. ID . To see that ( 16.31) is satisfied, we rewrite ( 13.11) as 



where we have used (1 - 



1 

2ik 



d_ 

dx 



(R r e)f 



— -Rr 

2ikdx n 



(6.4) 



(6.5) 



£ 2 )/ = (d/dx)i (equation (3.14)). Therefore, 

/ dzB\Rr(z,-oc;W;k)+Z(z,W)] = / dz -^-R r (z, -oo; W; k). 
J x-l 2ikJ x _ L dz 

The right-hand side vanishes because R r (x, — oo; W; k) is periodic in x. Thus R r + £ 
satisfies condition (16. 3p . It is not difficult to show, by using ( 13. ip and ( I2.15&P . that 
R r + £ is piecewise continuously differentiable in x and analytic in W on the real axis. 



7. Validity of the expansion 

Consistency between (16.31) and g = Ah requires that 

/ BA~ l g{z, W) dz = 0. (7.1) 

J x-L 

As shown in appendix G, equation ( 17. II) holds if g satisfies ( 15. 21) . This means that 
BA~ l g is in the domain of A" 1 if g is in the domain of A~ l . Recursively, it follows that 
(BA~ l ) n g belongs to the domain of A' 1 if g is in the domain of A -1 . In other words, 
A~ l {BA~ l ) n g makes sense for any n if A~ x g makes sense. 

New we can see that equations (13.111) . (13 . 1 2[) and (13 . 1 3[) are valid for periodic 
potentials if ( 13.31) is replaced by ( 14.31) . Since ( 13.12&P can be written as 

f n = 2Vr 1 (&4~ 1 ) n (l - e 2 )/ (n > 1), (7-2) 

it follows from the above argument that f n makes sense for any n if (1 — £ 2 )/ belongs 
to the domain of A^ 1 . From (I3.14p it is obvious that (1 — £ 2 )/ satisfies conditions ( 15.11) 
and (15. 2p . Therefore (1 — £ 2 )/ is in the domain of A^ 1 , and hence f and all f n make 
sense. It is easy to show that (^4 — 2ikB)~ 1 g makes sense if A~ x g makes sense. So the 
right-hand side of ( 13.131) makes sense, too. 

Thus, expansion (13. lip is justified for any N. The coefficients of the expansion f n 
given by ( 13.121) . as well as the remainder term p^, all make sense and are finite. 

The behavior of the remainder term as k — > is simple in the periodic case. 
From ( I3.13P and ( 14.1 ft we have lim^o pN I i}k) N+1 = A^Atn+i = fjv+i- Therefore, 
Pn = 0(k N+1 ) as k — > 0, since f n is finite for any n. Unlike in the non-periodic cases 
[2], no subtleties arise when the potential is periodic. Expansion (13.111) is asymptotic 
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for any N. Moreover, the infinite series ^^InO^)™ 7 ™ ^ s convergent for sufficiently small 
\k\. An explanation about the convergence of this series is given in appendix B. 

8. Expressions for f n 

Now let us calculate the coefficients f ,fi,f 2 , ... of (13.111) . First, f is given by ( 13.12ajl . 
Substituting (KIM into (03]), and using f*_ L e ±v{z,) dz' = L e ±v °, we obtain 



A-\l-i 2 )f = - . - — — V / dz' dzsxBh[V{z')-W\ — tanh 
v )J L smh(W-V ) J x -l J X . L 1 l dz 

— -i -X? r dz' sinh[y(z') - W] ( tanh - tanh 

3 inh(W-V ) A-l V 2 2 

-£> [L - L cosh(W - Vq) + L smh(W - V Q ) £(x, W)\ 



W-V(z) 



Lq sinl 

1 



L sinh(PV - Vq) 
sinh(iy - Vq) [1 " C ° SHW ~ Vo) + S[nHW ~ Vo) W)] 

W - Vn 

-tanh r— ^ + (8.1 



and hence 



W - Vq 

f = -tanh —5. (8.2) 

This is to be compared with (13. 16p . 

Next, we calculate = 2^4~ 1 i3(f + 0- We can see that 

9 f / Vn - W W - V(x) 
B (f + = \ sinh[PV - V»] (tanh + tanh ^ 

^ f • i rrr „, m / ,W-V ,Vo-V(x) 
= — |smh[V - V{x)\ (^tanh — + tanh - 1 } 

= sinh[V - V{x)\ 2co J w Vn - (8-3) 

Substituting this into (14. 3p . and multiplying by 2, we have 
1 

r 1 



L smh(W -V 



1 



cosh 



x V — \ 2 w-v 1 1 dzdz ' sinh[VV) - W] smh[V - V(z)}. (8.4) 



x— L<z<z'<x 



The integrals on the right-hand side can be expressed as 
II dz dz' sinh[VV) - W] sinh[V - V{z)\ 

J Jx-L<z<z'<x 

= \ (e- W+Vo i-+Tx- L + e W - Vo i+-T x „ L - e- w - v n++Tx- L ~ e W+V °[—]t_ L ) . (8.5) 
Substituting (I8.5P into (18.41) . and using definition (12.271) . we find 

1 r / e -w+Vo \ _ 2 

fl = 4L sinh(W-Vo)\ 1 cosh 2 ~ 1 J ( [ ~ + ^~ L ~ e " 2 
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a W-Vo 



cosh 



2 W-Vp 



1 -*"*[— M-l) ■ (8-6) 



2 

Using (12.221) and (I2.23p . we can reduce (I8.6P to the simple form 

f ^ 4LoCO I^ »+^-[- + ]^). (8.7) 

For n > 2, we can calculate f„ by successively applying C to fi as f n = L n ~ x f\. 
Let us write the operator B as 

b = \Y.^ v{x) ^ & = s w ( 1 + °m)- (8 ' 8) 

Substituting sinh[VV) - W) = ^J2a'=±i cr ' ea ' V{z ' )e ~ a ' W into (O, and using (EH}, we 
can write £ in the form 

C = 2A- 1 B = -L V V I^K-*,-*, (8.9) 
2L o CT=±1(J , =±1 

where we have defined the operators Z^/ and /C^o-/ as 

X a>al g(x,W) = if dzdz' e aV ^ +a ' v ^g(z,W), (8.10) 

</ Jx-L<z<z'<x 

-W 



sinh(Vb - W) 



/ dWX'X<?(z,W0- (8.11) 



The expression r 2 = Cr x with (18 .7p and (18. 9p reads 

= ^2 E E ' C0S h 2 W ) ( [+ ~^ " • (8 ' 12) 



Note that 



Iz-L 

= [vi, <T2)1-l + [<ri)*» + ki, (8-13) 

Therefore, 

J Jx-L<z<z'<x 

= [<7i, cr 2 , CT, </]*_ z + [cr 2 , £7, 01, ct']^_ l + [cr 2 , C7, O"', CTi]*_ L 

+ [<t,<ti,<T2,o^Z-l+ [^,o-i,o-',a 2 ]l_ L + [<T,a , } a u a 2 ]%-L- (8-14) 

This gives explicitly 

^++ ( [+-E-L - 1-+} x x -l) 

= [-+++)l-L + \+ + + -Tx-L ~ [+-++&-L ~ [++-+)l-L 

= 1P 3 M-P[+-+]^_ L , (8.15a) 
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*-([+-]*-£- [-+K-J 

= -[+ ]Li - [ +]%-L + + [-+—]'*-L 

= -lM 3 P + M[-+-]:_ L , (8.156) 

If- ( [+-M-L - \-+]l- L ) = l-+-+)t-L + [+-+-)t-L ~ 2 [+—+M-L 

= 2Q-M[+-+}:_ L} (8.15c) 

X-+ ( [H — IS-zi - = -[-+-+\Z-l - [+-+-1^ + 2 

= -2Q + P[-+-]Z_ L , (8.15d) 

where we have used (I2.20p and defined 

Q = [-+-+]Z-l+[+-+-]1-l- (8-16) 

This Q is independent of x, as can be seen by differentiating the right-hand side with 
respect to x and using ( 12.24]) . The expression in large parentheses in f)8.12p can be 
calculated as 

I a / e -(a+a')Vo / e -[{a/2)+a'}(W-Vo) 

K ^-°' C o S h 2 ^ = sinh(^ - V ) { cosh 3 ^ 1 ' ' (8 ' 17) 
Substituting flSTTS) and flSTTj) into (jSJ) , we obtain 



4L n cosh 3 ^> 



The coefficient r n for n > 3 can be obtained in the same way. The expression for 
general n is 



"n 2 n+l L n ^ ( ' ' ' 

\ 



^' — -^ cosh 2 W^fc 



x ^..C, • " -Z^Z^ ( [+-] X X -L ~ , (8-19) 

where the sum is over <jj = ±1 and o\ = ±1 for 1 < i < n — 1. 

9. The expansion of # 

It was shown in {2j that the expansion of S r takes the form 

S r (x, k) — - — a (x) + \kax{x) + (iA;) 2 a 2 (x) + (iA;) 3 a 3 (x) + • • • , (9.1) 



where 



4 W-4—oc 



^Jim e- w+v ^(f -l), (9.2a) 
a n = \ lim e~ w+v{x) f n (n > 1) (9.26) 

4 
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(see equations (4.4) of 0). Substituting (Q, dHUJ) and (KTTgj) into (0 gives 

ao(x) = -ie v ^- Vo , (9.3a) 

^)-^ w {^i+^--5z;(t + «)}- (9 - 3c) 

When the potential is periodic, the functions SV and Sj are simply related by 

3 r (x,k) = 3i(x,-k). (9.4) 
(See appendix D. This is a special feature of the periodic case.) Therefore, 

S(x, k)-l = s (x) + (ik) 2 s 2 (x) + (i£0 4 s 4 (x) + • • • , (9.5) 

where 

s n ( x ) — 2a„(a;) (n even), s n (;r) = (n odd). (9.6) 
The first two nonvanishing coefficients are 

So (x) = -e y ^- y °, (9.7a) 

10. The expansion of the Green function 

The small- expansion of the Green function is obtained by substituting ( 19 .51) into (12.131) . 
The expansion has the form 

G a (x, y; k) = + a + iAtfi + {ik) 2 g 2 + ■■■, (10.1) 
and the coefficients can be expressed in terms of Sq, si, si, . . . as 

= , \ , (10.2a) 
2A/s (x)s (y) 

9o(x, y) = qi(x, y)g-i(x, y), (10.26) 

»(*,») = l{ [qi{ *,y)r-^-^} M , (10.20 



9z(x,y) = \ hqi(x,y)] 3 + q 3 (x,y) - -qi(x,y) 



S2(x) s 2 (y) 



g^{x,y), (W.2d) 



6 2 [s (x) s (l/). 

where we have defined 

q n (x,y) = - s n _i(z) dz. (10.3) 

•A/ 

Substituting (19. 7p into (110. 2p . and also using q±(x,y) = e~ v °[+]y, we obtain 

e -[V(x)+V(i0]/2 

^_i(x,3/) = e^ , (10.4a) 
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[V(x)+V(y)}/2 

9o(x,y) = [+]* (10.46) 

e -[V(x)+V(y)]/2 f 

s^y) = — — ( [+-+]*-£ + 

+ ^oe-^([+B 2 -^(^ + g)}, (10.4c) 



e- 3Vo , r .v,xs 



g 2 ( x ,y) = e- v °{+r ygi (x,y)-\ — (i+];y+ I s 2 (z)dz\g^(x,y). (10 Ad) 



The series f llO.ip is convergent for sufficiently small \k\ (see the comment at the end of 
section 7 and appendix B). 

11. Example 

As an example, let us consider the periodic square potential (22] 

f (0 < x < a) 
V(x) = { / V(x + L) = V(x). 11.1 

I C (a < x < L) 

For < x < a, we have the exact expressions of a and (3 as [17J 

a(x, x-L;k) = —t— e - ikL (l - A 2 e 2ikb ) , (11.2a) 
1 — A 2 

f3(x,x- L; k) = A e 2^ e -ifcL / e 2ifc6 _ 1 \ (\\.2b) 
1 — A 2 

where 

A = — tanh(C/2), b = L - a. (11.3) 

From (lll.2p we obtain R r , Ri and S by using ( ID. 71) and ( 1D.9I) of appendix D. The 
expressions for a < x < L have similar forms. The exact expression of the Green function 
is obtained by substituting ( ID. 91) into ( 12.1 3D . or, more simply, by using equations (3.6) 



of [20]. The resulting expression of G$ for 0<y<x<ais 

_ [2Ae 2ikx e- ik{L - b) sin kb - K(k)] [2Ae- 2 ' lky e lk{ ~ L -^ sin kb - K(k) i 

Gs(X ' m ] = 2ik^-v) {4A 2 sin 2 kb-[K(k)f) ' (1L4) 
where 

K(k) = smkL - A 2 sin[k(L - 2b)} - (1 - A 2 ) {l - [F(A;)] 2 } 1/2 , (11.5) 

Y(k) = 1 {cos kL - A 2 cos[k(L - 2b)] } (11.6) 

(see (15.81) for the definition of Y). 

It is characteristic of periodic systems that the energy spectrum has a band 
structure. Let us suppose that A; is a real number. The energy lies in a band if 
|Y(/c)] 2 < 1, and in a gap if [Y(fc)] 2 > 1 (see appendix D). The band structure can 
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Figure 1. Graphs of (a) Y(k) 7 (b) Re Gs(x, y; k) and (c) ImGs(x,y;k) for the 
potential (jll.lj) . plotted as functions of real k, with C = 1, L = 1, a = 0.6, ir = 0.4 
and y = 0.1. In (a), the shaded and unshaded areas represent the bands (\Y\ < 1) and 
the gaps (|F| > 1), respectively. In (b) and (c), the solid lines show the exact value of 
Gs (equation (|11.4jl ). The imaginary part of Gs is identically zero in the gaps. The 
graphs on the right are close-up plots, where the dashed lines show the result of the 
low-energy expansion up to order k 2 . The dashed curves almost overlap with the solid 
curves in the lowest band except for the imaginary part near the band edge. 



be seen in the graph of Gs{k) (figure 1). The low-energy expansion can be used for 
calculating the Green function in the lowest band. 

In this example, M = a + be~ c and P = a + be c , as can be easily seen. Hence, 



^J{a + be- c ){a + be c ), 



3 Vb 



a + be c 



;n.7) 



a + be~ c 

We assume that < x < a. The integrals appearing in the expansion of R r or 5* can be 
calculated, for example, as 

[+-F-L = [+-1S + [+-]% + 1+-V-L + [+V- L H- b + [+V-lHo + 1+toHo 
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= - [x 2 + b 2 + (a - x) 2 ] + e~ c (a - x)b + (a - x)x + e c bx 

= ^(a 2 + b 2 ) + e~ c b(a- x) +e c bx. (11.8) 
Changing the sign of C gives 

[-+M-L = \(a 2 + b2 ) + e °b(a - x) + e~ c bx. (11.9) 
In a similar way, 

[+-+]*, = -(a 2 + 36 2 ) + 2&x(x-a)sinhC + -(3a 2 + 6 2 )e c , (11.10a) 
6 6 

[-+-]* = ^(a 2 + 36 2 ) - 26x(x - a) sinhC + ^(3a 2 + 6 2 )e" c , (11.105) 
6 6 

= ^(« 4 + 6a2 b 2 + b 4 ) + ^[b(3a 2 + b 2 )x - 6abx 2 + 4bx 3 } sinhC 

+ y(a 2 + 6 2 )e" c , (11.11a) 
[-+-+]*-£ = ^(« 4 + 6fl2&2 + &4 ) - \HZcl 2 + b 2 )x - 6abx 2 + 4bx 3 } sinhC 

+ ±(a 2 + b 2 )e c . (11.116) 
o 

Substituting (lll.lip into (JEHU), we nave 

Q = —(a* + Qa 2 b 2 + 6 4 ) + —(a 2 + b 2 ) coshC. (11.12) 
12 3 

For < y < x < a, it is obvious that [+]* = x — y. The low-energy expansion of 

Gs(x,y;k) for < y < x < a is obtained by substituting the above expressions into 

( 110. 4p (and (19. 7} ). The result of the expansion up to order k 2 is plotted in figures 1(b) 

and (c) (the graphs on the right). It can be seen that this expansion gives a very good 

approximation in the lowest band except near the band edge. 



12. Conclusion 

The low-energy expansion formula for the reflection coefficient introduced in [TJ and [2] 
for non-periodic potentials (equations (13. lip . (I3.12p and ( 13 . 1 3[) of the present paper) is 
applicable to periodic potentials as well, if only we replace (I3.3P by (14. 3p . In the periodic 
case, the general expression for the nth-order coefficient of the expansion is obtained 
as ( I8.19p . The expansion of the Green function takes the form of (llO.ip . with the first 
four coefficients given by (I10.4[) . The method of this paper can be extended to more 
general potentials which are asymptotically periodic at spatial infinity, and this will be 
discussed elsewhere. 



Low-energy expansion formula for periodic potentials 24 
Appendix A. Derivation of ( 12.101) 

As noted in section 2, the functions a(x,x';±k) + (3(x,x';±k) are solutions of the 
Schrodinger equation (jl.3p . If / is identically zero (i.e. V = constant), then obviously 
a(x,x'; ±k) = e^ lk ^ x ~ x "> and (3(x,x'; ±k) = 0. Hence, 

a{x,x';±k) + f3{x,x';±k) =e Tik{x ~ x>) if / = 0. (A.l) 

Note that definition (12.31) is not restricted to x > x'. Since U(x', x; k) is the inverse of 
the matrix U (x, x'; k), we have 

a(x',x;±k) = a(x,x',^fk), f3(x', x; ±k) = — (3(x, x'; ±k). (A. 2) 

Let us replace f(x) by f XljX2 (x). If we define 

Xi(x) = a(x, Xi, k) + p(x, x\\ k), (A. 3) 

then Xi( x ) is a solution of the Schrodinger equation (|2.9p . Since f XltX2 (x) — for x < xi, 
it is obvious from ( 1A.1I) that 

Xi(x) = e ~' lk{x ~ Xl) for x < x 1 . (A.4) 

From (12. 5p we have 

a(x, Xi,k) = a(x, x 2 ; k)a(x 2 , X\\ k) + f3(x, x 2 ; —k)/3(x 2 , X\\ k), 

(3(x, xi,k) = (3(x, x 2 \ k)a(x 2 , x\\ k) + a(x, x 2 \ —k)f3(x 2 , x\\ k). (A. 5) 

Using this, we can write 

Xi{x) = a(x 2 ,xi; k)[a(x,x 2 ; k) + /3(x,x 2 ; k)) 

+ f3(x 2 ,Xi, k)[a(x,x 2 ; -k) + f3(x,x 2 ; —k)]. (A. 6) 

Since f Xl , X2 (x) = for x > x 2 , from (1A.1I) and (1A.6I) we find 

Xi{x) = a{x 2 ,x 1 ;k)e- [k{x - X2) +P{x 2 , Xl ;k)e ik{x - X2) for x > x 2 . (A.7) 

If we divide the right-hand sides of (1A.4|) and (IA.7[) by a(z2,^i; fc), they coincide with 
the right-hand sides of ( I2.10al) . Hence, we can see that ipi(x) = xi(x)/a(x 2 , x\\ k) is the 
solution of the Schrodinger equation (12. 9ft satisfying (12. 10 aft . 
In a similar way, we define 

X%{x) = oc(x, x 2 ; -k) + /3(x, x 2 \ -k). (A. 8) 

By using 

a(x, x 2 ; —k) = a(x, x\\ —k)a(x 2 , x\]k) — (3(x, x 2 ; —k)/3(x 2 , x\, —k), 
(3(x, x 2 ; —k) = (3(x, x\\ —k)a(x 2 , xi,k) — a(x, x 2 ; —k)(3(x 2 , x\\ —k), (A. 9) 
we can rewrite X2 as 

X2(x) = a(x 2) x\] k)[a(x,xx; -k) + /3(x,xi, —k)] 

— (3(x 2 , xi, — k)[a(x, Xi, k) + (3(x, x\\ k)). (A. 10) 

It is obvious that ip 2 (x) = X2(x)/a(x 2 ,Xi, k) is the solution of the Schrodinger equation 
satisfying (12101 . 
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Appendix B. Analyticity of ot(k), (3(k), and the reflection coefficients 

We define the functions a^ 1 and by 

a+(x,x';k) p-(x,x';k) \ = ±( 1 1 \ f a ( k ) P(~ k ) \ I 1 ~ l 
/3 + {x, x 1 ] k) a-(x,x';k) J = 2 I -1 1 J I f3(k) a(-k) ill 1 

if a(k) + a{-k) + p{k) + p(-k) ~a{k) + a(-k) - 0{k) + p{-k) 
' 2 I -a(k) + a(-k) + (5{k) - (5{-k) a(k) + a(-k) - f3(k) - p(-k) 
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(B.l) 



where a{k) and (3(k) stand for a(x, x'; k) and fi(x,x'; k). From (12. II) . it follows that 



d { a + (x,x';k) (5 (x,x';k) 
dx \ /3 + (x,x';k) a~{x,x';k) 

Using f ll.2p . we can rewrite flB.2|) as 

i v W 2 a + (x,x';k) 



f(x) ik 
ik -f(x) 



a + (x,x';k) (3 (x,x';k) 
f3 + (x,x';k) a~(x,x';k) 



(B.2) 



d 




dx 


e 


d 




dx 


e 



ike v ^ 



- v M' 2 P + (x,a/;k) 



v ^ 2 P + (x, x'; k)} = ike~ v ^ \e v ^ 2 a + (x, x'; k) 



(B-3) 



The initial conditions for these equations are a + (x',x';k) = 1, j3 + (x',x';k) = 0. 
Equations ( IB.3j) with these initial conditions are equivalent to the integral equations 



i V ^ 2 a + {x, x'- k) = e v(x 'V 2 + ik / e v « \e- v ^ 2 p+(z, x'; k) 



dz, 



dz. 



e- v W 2 p+(x,x';k)=ik / e~ v ^ e v ^ 2 a+(z, x'\ k) 

Jx' 

These equations can be solved by iteration. The solution is formally expressed as 

oo 



(B.4) 



m=0 

oo 



(B.5) 



m=0 

where l% m and I^m+x are defined recursively by 

ILM= r* VW lL-iMA*, tin+iM= f e- y ^lt m {z,x')dz, (B.6) 

Jx' Jx' 

with Iq = 1. It is easy to see that I^(x,x') = [—]%>, a/) = [ — \-}*/, I^ix^x') = 

[ — I — Yx'i anc ^ so on - Namely, 

Itjx.x 1 ) = [-+-+..._+]», lL + i(^x') = [_+-+...+_]». (B.7) 

In the first equation of (IB. 71) . there are 2m alternating plus and minus signs in the 
expression on the right-hand side. In the second equation, there are 2m + 1 alternating 
signs. 
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We assume that both x and x' are finite. As we are assuming that V(x) is finite 
for any finite x, there exists a constant number C such that 

e _ y(z) < ^ eV{z) <(J for x r < z < x ( R g) 

From definition ( I2.19p . we can easily see that 

I 2m (x, x ) < ^2 m )i ' hm+i( x i x ) < (2m + 1)!' 

Therefore, 

E |w 2m 4n(^')| < E i fc i 2m (2^)! = cosh (^i) > 

m=0 m=0 

E K ifc ) 2m+1/ 2 + m+ i(^^)| < E l fc | 2m+1 f2 ^ + 1) , = *nh(C\k\) . (B.10) 

m=0 m=0 ^ 

Thus, the infinite series on the right-hand sides of ( 1B.5|) are absolutely convergent for 
any fc. This means that fe) and (3 + (x,x; k) are entire functions of k. 

In the same way, we have 

oo 

a-(x,x'; k) = e ™- v VW 2 E(^) 2m/ 2m(^^')> 

m=0 

oo 

r(x ) x^) = e-^) + ^/ 2 E(^) 2m+1 W 1 (^^') ) (B.H) 

m=0 

where 

„/\ r i i , la; t— /„ „,A 



J 2 - m (x, x') = [+-+-•■ ■ +-]*, i^ n+ i(a;,x') = [+-+ +]£. (B.12) 

Equations ( IB. Ill) and (IB. 121) are obtained from (1B.5I) and ( IB .71) by changing the sign of 
the potential V. Just like a + and $ + , the functions a~ and $~ are entire functions of 
k. Inverting (IB.ip . we can express a(±k) and /3(±k) in terms of (X s1 and as 

a{x,x';k) 0(x, x'\ -k) \ _ 1 / a + + or - /3 + - f}~ a + - oC - /3 + + f}~ 
fi(x, x'; k) a(x,x';-k) J ~ 2 1 a + - or + f3 + - j3~ a + + or + /3 + + 0' 

'(B.13) 

Obviously at(a;, a;'; fc) and j3{x,x'\ k) are entire functions of k as long as both x and x' 
are finite. We obtain (15. lip by substituting ( 1B.5[) and (IB. lip into (IB. 13[) . 

It is shown in appendix D that the reflection coefficient i? r (x, — oo; k) (\mk > 0) 
for periodic potentials can be expressed in terms of a(k) = a(x, x — L;k) and 
(3{k) = f3(x, x — L;k) as 

*<*--^ t) = [i/A(*M-U) (R14) 



with A(fc) = Y(k)-i^l - [Y(k)} 2 and Y(k) = [a(k)+a(-k)]/2. Since a(Jfe) and f3(k) are 
entire functions, R r (x, — oo; fc) may possibly have singularities only where 1/X(k) = a(k) 
or Y(k) = ±1. 
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It can be shown that \r(x2,x\] k)\ < 1 for Imk > and X\ < X2 (see section XIII 
of p2]). This implies that |a(A;)| > 1 for Imk > 0. Since \X(k)\ > 1 for Imfc > 0, the 
equality 1/X(k) = a(k) can never hold in Imk > 0. If Y(k) = ±1, then |A(/c)| = 1. This 
cannot hold for ImA; > 0, either. (If |A(A;)| = 1, there exist two independent solutions 
of the Schrodinger equation which remain finite as x — > ±00. This can happen for 
Imk = 0, but not for Imk > 0.) Thus, neither 1/X(k) = a(k) nor Y(k) = ±1 can hold 
in Imk > 0, and so R r (x, —00; k) is analytic in the upper half plane (Imk > 0). Note 
that R r (x, —00; k) has branch point singularities on the real axis. These branch points 
correspond to the edges of the bands. 

The power series expansion of R r (x, —00; k) in terms of k is convergent for 
sufficiently small \k\ if the right-hand side of (IB. 141) is analytic in some neighborhood of 
the origin. Since 

[1/X(k)] - a(k) = 2ikL cosh 2 ~ V ° + (k 3 ), 

l-[Y{k)} 2 = k 2 L 2 + O(k 4 ), (B.15) 

there exists a number C such that [l/\(k)]—a(k) 7^ and 1 — |F(/c)| 2 7^ for < \k\ < C. 
Although [l/A(/c)] — a(k) = at k = 0, this zero is canceled by the zero of (3{k) in the 
numerator (see equation (15.161) ). and so the right-hand side of (IB. 140 does not have a 
pole at k — 0. The origin is not a branch point, either, since y/l - \Y(k)\ 2 ~ kL for 
k ~ 0. Therefore R r (x, —00; k) is analytic, and the power series expansion in terms of 
k is convergent, in \k\ < C. 

Since R r (x, —00; W; 0) = — tanh ~ 7^ 00 (see (18.21) ). we know that R r does not 
have a pole at the origin. As R r (x, — 00; k) is analytic in a neighborhood of k = 0, 
it is obvious from definition (12.1561) that R r (x, — 00; W; k), too, is analytic in some 
neighborhood of the origin. Thus, the expansion of R r (x, —00; W; k) in terms of k is 
convergent for sufficiently small \k\. 



Appendix C. Relation with the Weyl-Titchmarsh m-function 

A natural way to define the transmission and reflection coefficients (for finite intervals) 
for the Schrodinger equation is to consider the truncated potential 

{0 (x < X\) 

V s (x) (x 1 <x<x 2 ) (C.l) 
(X2 < x). 

This Vjf is different from f% ix + f' xlX2 (which appears in (12.91) ). since the latter 
contains delta functions at x = X\ and x = x 2 (unless f(xi) = or f(x 2 ) = 0). The 
Schrodinger equation 

- S + o> = ^ < c ' 2 > 
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has two independent solutions of the forms 

, / \ _ I T S (x2,xf, k) e" ifc (— (a; < Xl ) 

e ik ^ + Rf{x 2 , x i; k) e~ ik ^ {x < x x ) 
' 2 '"' i ^ r s (x 2 ,x i; A;)e ifc( ^ 2) (x>x 2 ). { ' } 

This defines the transmission and reflection coefficients r s , Rf and Rf for the 
Schrodinger equation. Here we assume that Im k > 0. When Vs belongs to the so-called 
limit point case [7], we can let X\ — > — oo or x 2 — > +oo to obtain the reflection coefficients 
for semi-infinite intervals. As shown in [2], the reflection coefficients Rf(x, — oo; k) and 
-Rf (oo, x; k) are related to R r (x, — oo; k) and Ri(oo, x; k) by 

R r (x,—oo;k) Rf(x,—oo;k) f(x) 

1 + Rr(x, — oo; fc) 1 + Rf(x, — oo; fc) 2ifc 

Ri(oo,x;k) Rf(oo,x;k) f(x) 

+ ^777- 



1 + Ri (oo, x; k) 1 + Rf(oo,x;k) 2ik 

Let 771 and r/ 2 be two independent solutions of the Schrodinger equation (11.31) 
satisfying the conditions 

rfr(x ) = l, r][(x o ) = 0, r] 2 (x ) = 0, t/^(x ) = 1, (C.6) 

with some fixed xo- And let x± denote the solutions of (11.31) satisfying x+(°°) — and 
X-(-oo) = (and which are are square integrable near +00 and —00, respectively). 
These functions can be expressed as linear combinations of 771 and 772- We write 

x± (x) = %(x) ±m ± r] 2 (x). (C.7) 

This is the definition of the Weyl-Titchmarsh functions m + (xo, k) and itl_(xq, k), which 
are regarded as functions of xq and k. 

Let us suppose that Vs(x) = for x > x . Then, 

771 (x) = cos[A;(x — xq)], f] 2 (x) = — sin[A;(x — xq)] for x > xo, (C.8) 

K 

and so 

X-(x) = cos[/e(x — xq)\ sin[fc(i — xq)\ 

K 

= ^-[k + im_(x , k)} e ik(x - x °1 + J-[Jfe - im_(x , k))} e~' lk ^ Xo) for x > x . (C.9) 
Comparing this with (IC.3I) . we find [231 E 



dS/ u \ k + im_(x ,k) Rf(x ,-oo;k) 

Rr(x ,-oo;k) = - : — , m_(x ,k) = \k-2\k — . (C.10) 

k — vm-{Xo, k) 1 + Rr(xo, — 00; k) 

In the same way, 

oS/ n k + im + (x ,k) Rf (00, x ; k) 

R?(oo,x ;k) = : — , m+(x ,k) = lk - 2ik s — . C.ll 

k — im + (xo, K) 1 + Rf (00, x ; k) 
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From (O, (1U1Q|) and (lUTTj) . we obtain 

5 r (x, fc) = ^[m_(x, fc) + /(*)] + ^, fc) = ^[m+(x, k) - f{k)\ + ± (C.12) 
and 

5(a;, fc) = ^[m+(z, fc) + m_(x, fc)] + 1. (C.13) 

Appendix D. Exact expressions of R r , R[ and S 

The function a(x 2 ,xi; k) satisfies the differential equations 
d 



-a(x2, x\] k) = —ika(x2, X\\ k) + f(x2)/3(x2, x\\ k), 



dx 2 
d 

- — a(x 2 , xi, k) = ika(x 2 , Xi, k) - f(x 1 )/3(x 2 , X\\ -k). (D.l) 

OX i 

(The first equation of (ID. 1|) is a component of (12. ip . The second equation is obtained 
by using (|A.2[) of appendix A.) Hence, if / is periodic, 
d 

—a(x, x — L;k) — f(x) [f3(x, x — L; k) — (3(x, x — L; —k)) , (D.2) 

and so 

d 

— \a(x, x — L; k) + a(x, x — L; —k)] = 0. (D.3) 
ox 

Thus, Y{k) defined by (15. 8p is independent of x. We define 

Z(k) = {l-[Y(k)f} 1/2 (D.4) 
to write (15. 7)) as 

A (A;) = Y(k) - iZ(k), 1/ A(fc) = Y(k) + iZ(k). (D.5) 

The branch of the square root in (1D.4|) is chosen such that \\{k)\ > 1 for Imk > 0. 
With this choice of the square root, Z(k) is an odd function of k whereas Y(k) is an 
even function (see ( 15. 13[) ) . So, from (ID .5|) we find 1/X(k) = X(—k). 

Hereafter we write a(k) and (3(k) in place of a(x,x — L;k) and (3(x,x — L;k), 
respectively. Substituting the elements of (I5.9P into (12.76ft gives 

-/9(fc)(l- 7 ») 



R r (x, x — nL; k) 



[l/\{k)\ - a(k) - -f n [\(k) - a(k)} 



Rl (x + nL, x- k) = Rl (x, x - nL; k) = 1-^1-^ _ ^ , (D.6) 

where 7 = [A(/c)]~ 2 . Since I7I < 1 for Imk > 0, we can take the limit n — > 00 of these 
expressions and obtain 

-f3(k) 2/3(k) 



R r (x, —00; k) 



[\/\{k)]-a{k) a(k)-a(-k)-2iZ(k) 



ff(-fc) _ -2P{-k) 
X] } ~ [1/X(k)} - a(k) ~ a(k) - a(-k) - 2iZ(k) ' { } 
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Substituting into (12.111) . and also using ( 12. 4p . we have 

2/3(k) la(k)-a(-k) + 2/3(k) + 2iZ(k) 



Si(x, k) 



a(k) - a(-Jfe) + 2(3{k) - 2iZ(k) 2 a{k) - a(-k) + (3{k) - p(-k) ' 

-2p(-k) 1 a(k) - a(-Jfe) - 2p(-k) + 2\Z{k) 



a(k)-a(-k)-2p(-k)-2iZ(k) 2 a(k) - a(-k) + p(k) - p(-k) 

(D.8) 

Since Z(—k) = —Z(k), it is obvious that Si(x,k) = S r (x,—k). The expression of S is 
obtained from (ID. 81) as 

S(x,k) = l + — . ~r~^~ryr — — . (D.9) 

v ; a(A) - a(-k) + p{k) - P(-k) 

By definition, the above expressions hold for Im k = as well. 

When k is a real number, [a(k)\* = a(—k), [fi(k)]* = (3{—k), and Y(k) = Rea(k). If 
[F(fc)] 2 > 1, the two eigenvalues of U(x,x — L; k) are real, and \X(k)\ > 1. If [^(A:)] 2 < 1, 
then |A(/c)| = 1. This means that the energy lies is in a band if [F(/c)] 2 < 1, and in a 
gap if |Y(A;)] 2 > 1 (see, for example, [7]). 

Appendix E. Details of the calculation leading to ( 15.201) and ( 15.211) 

To simplify the expressions, we write 
_ V{z) - V(x) _ Vjz) + V(x) _ V(x) - V _ W- V(x) 

s = 2 , * = - 2 , e = — - — , <p = — - — , (E.i) 

and 

b = \ ( [-+} X x -l ~ ■ (E-2) 
Then, (loTTTj) and (15TT8|) read 

a(x, x - nL- ±k) = (cosh 2 6) X ±n - (sinh 2 6) A T " + 0(A; 2 ), 

/3(x, x - nL- ±k) = ( T sinh cosh - iJfe-^- J (A n - A~ n ) + 0(k 2 ), (E.3) 
a(x - nL, z; ±k) = cosh 5 T j (e w Hr' + e"" [+]r™ L ) +0(A; 2 ), 

/9(x - nL, 2; ±Jfe) = sinh 5 ± y (e" [-]r™ L - [+} x z ~ nL ) + 0{k 2 ). (E.4) 

Substituting (jE3|) and (lE4j) into (loTT^]) yields 

a(x, 2; ±Jfe) = cosh cosh(0 + 5) X ±n - sinh sinh(fl + 5) X Tn 
T j coshtf (e^ 9 [-]r nL + e e "" {+f z ' nL ) A ±n 
qp y sinhfl (e^- 9 [-]f~ nL - e e ~^ [+]r™ L ) A Tn 

-iA;-^-sinh5(A n - A~ n ) +0{k 2 ), 
2L 
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0(x, z; ±k) = - sinh 9 cosh(# + 5) A ±n + cosh 9 sinh(fl + 5) A T " 
± j sinh# (e^ 9 [-} x z ~ nL + e 9 -" \+]T nL ) A ±n 

± y coshfl (e"~ e {-f z ~ nL ~ e'"" [+} x z ~ nL ) A Tn 

-iJfe— cosh<f(A"-A- n ) +0(A; 2 ). (E.5) 
2L 

Equations (IE.5I) can be written in the form 

a(x, z; k) f3(x,z;—k) \ I cosh# — sinh# \ / A(k) B(—k) 
(3(x,z;k) a(x,z;—k) J \ — sinh0 cosh# J I B(k) A(—k) 

cm I M !i ] . (e.6) 



where 



cosh 5 sinh 5 



A{±k) = cosh(fl + <T)A ±n =f - (e^ -0 Hr™ L + ^ i+lT^) A± ™> 
B(±Jfe) = sinh(# + 5)X Tn ± - {e u ~ 9 [-\ x ~ nL ~ e 9 "" [+]'~ nL ) ^ 
C(*0 = (A n - X- n ) . (E.7) 



2L 
Therefore, 

a(x,z;W;k) j3(x,z;W; —k) 
(3(x,z;W;k) a(x, z; W; — k) 



cosh<£> — sinh<£> \ / cosh# — sinh^ \ / A{k) B(—k) 
— sinh Lp cosh<y9 / I — sinh^ cosh# / I B(k) A(—k) 



-C{k) 



cosh cp — sinh cp \ I sinh 5 cosh 5 
— sinh cp cosh cp J \ cosh 5 sinh 5 



cosh(fl + ^) -sinh(fl + ^) \ ( A(k) B(-k) 
-smh(9 + ip) cosh(# + y?) J \ B(k) A(-k) 

V cosn(o — ip) smh(o — cp) I 

Comparing (IE. 61) with (1E.8I) . we can see that a and /3 are obtained from a and /3, 
respectively, by the replacements 

9 ->■ 9 + cp, 5 ->■ 8 — cp, u -> u + cp. (E.9) 

Namely, 

a(x, z; W; ±k) = cosh(fl + cp) cosh(fl + 5) A ±n - sinh(fl + cp) sinh(fl + 5) A T " 

iA- 

=F - cosh(fl + cp) (e"- 9 {-] x z - nL + e 9 "" \+} x ~ nL ) \ ±n 
\k 

=F - sinh(fl + y>) (e"" e [-]r nL - e^" [+]r nL ) >^ n 
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- ik— sinh(5 - <p) (A" - A" n ) + Oik 2 ), 

Pix, z; W; ±k) = - sinh(fl + <p) cosh(fl + 5) X ±n + cosh(fl + <p) sinh(fl + 5) X Tn 
± ^ sinh(# + if) (e w ~ 6 [-\ x ~ nL + e'"" [+]*~ nL ) A± ™ 



2 

±^cosh(# + ^) (e w ^[- 



x—nL g6>— u> r_|_~| — nL\ \--n 



: ) A Tr ' 



ik— cosh(5 - ip) (A n - A" n ) + 0(k 2 ). 
2L 



(E.10) 



Note that the replacements (1E.9|) are equivalent to V(x) —> W. This is an expected 
result if we consider the meaning of the variable W (see |21j). Equations ( I5.20p are 
obtained by substituting ( lETTj) and ( lR2j) into dETTOj) . 
From ( IE. 101) . we obtain 

a(x, 2; W; k) + /3(x, 2; W: -k) = e e+s (cosh(fl + <p) X n - sinh(fl + <p) X~ n ) 

- ike 1 "- (cosh(0 + <p)X n + sinh(# + tp)X~ n ) [-} x z ' nL 

- ' lk wr ^ ( A " " A_n ) + ( fc2 )' 

a(x, z; W; k) - fi{x, z; W; -k) = e~ {9+5) (cosh(0 + <p) X n + sinh(# + if) A" n ) 

- ike e ~ u (cosh(# + tp)X n - sinh(# + ^)A" n ) [+} x ~ nL 

+ ik^-^- & (X n -X- n ) +0(k 2 ). (E.ll) 

Substituting (lE.lj) . and using definitions (I5.22p . we can write ( IE. lip as 
a(x, z- W- k) + j3(x, z- W; -k) = e [v(z) - Vo]/2 cosh (1 + c ol n ) X n 

x ( I - ifc 1 " ^ e y " [-]r nL + kbA + 0(A; 2 ), 



1 + c 7" 

a(x, 2; VT; k) - p(x, z; W; -k) = e -[^)-vb]/2 cogh w-^ ^ _ A 

. , 1 + c 7' 



x 1 - ifc- 



1 - Co7" 

(E.12) 

where the quantities b\ and 62 (which come from the terms involving h) are independent 
of z. Equations ( 15.211) are easily obtained from ( IE. 121) . 



Appendix F. Verification of AA 1 g = g and A 1 Ah = h 

Let us verify AA~ x g = g and A~ x Ah = h using expression (14.3)) . or equivalently ( 15. 33j) . 

We can see that 
d f x . 



, ix _ Q Wi ix 



r z )g(z,W)dz 
e w [-r x _ L )g(x-L,W) + (e v ^- w 



) / g(z,W)dz 



x—L 



2L smh(W-V )g(x,W), 



(F.l) 
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where we have used (I2.24p . (15. ip and (15. 2p . Since the last line of (IF.ip vanishes for 
W = Vq, we can insert the operator T> in front of the integral in the first line without 
changing the result. Namely, 

AV [ (e- w [+} x z -e w [-} x z ) g(z,W)dz = 2L smh(W -V )g(x,W). (F.2) 

Jx-L 

Dividing both sides by 2L sinh(W — Vo), we obtain AA~ l g = g. 
In a similar way we have, by integrating by parts, 

£ L (e- w [+)*-e w [-]*)^h(z,W)dz 

= 2L sinh(jy - V )h(x,W) -2 / sinh[W - V(z)]h(z, W) dz. (F.3) 

Jx-L 



We apply the operator V to both sides of (1F.3|) . The first term on the right-hand side 



does not change. For the second term, we can use (12.281) to write 

rx rW i-x o 

V / sinh[W -V(z)]h(z,W)dz = / dW dz — - sinh[W - V(z)]h(z, W) 

Jx-L Jvo Jx-L 



V 

W rx 

dW \ dzBh(z,W). (F.4) 

V Jx-L 



This vanishes on account of (16. 3p . Therefore, 

v! (e- w [+] x z -e w [-} x z )Ah(z,W)dz = 2L smh(W-V )h(x,W). (F.5) 

Jx-L 

Dividing by 2L Q sinh(PV — Vo) gives A _1 Ah = h. 

Appendix G. Verification of ( 17. lj) 

We can write the left-hand side of (17.11) as 

/ BA~ 1 g(z,W)dz= } T Q f T/ ( . * T _, / dzsinh[PV - V(z)} 
Jx-L 2L dW \ smh(PV - Vo) J X _ L 1 n 

x vj' M {e- w [+)i -e W [-)i)9V, WO}- (G.l) 

Substituting 

sinh[W - V{z)\ = 1 A ( e -^[+]; - e w [-Q , (G.2) 
and integrating by parts, we obtain 
f X dzsmh[W -V(z)} f Z dz' (e~ w [+] z z ,-e w [-]*,) g(z',W) 

J x—L J z—L 

= -l(e- w P-e w M) f dz'{e- w [+Y z ,-e w [-]*,) g{z\W) 

1 Jx-L 

+~ £ dz (e- w [+} x z - e w [-}:) (e~ w P - e w M) g(z, W) 

- f" dz(e- w [+] x z -e w [-] x )smh[V(x) -W] f dz'g(z',W) 

J x—L J z—L 

= 0, (G.3) 
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where we have used (15.21) . In a similar way, using e W P — e w M = 2L sinh(Vo — W) 
and e~ v °P - e v °M = 0, we have 

f dz sinh[ty - V(z)} f dz' (e" y ° [+]*, - e v ° [-]%) g(z', V ) 

J x—L J z—L 

= L smh(W-V ) [ X dz' (e- Vo [+r z , - e Vo [-y z ,) g(z',V ). (G.4) 

Jx-L 

Equations ( 1G.3I) and ( 1G.4I) give 

smh(W - V ) [ l ^ Shim ~ V(Z)] V f l^' ( e ~ W[+]l ' ~ eW[ ~ ] ^ 9{A W) 

= lJ X dz'(e- y °[+F z ,-e v °[-F z ,)g(z',V ). (G.5) 

Jx-L 



Since this expression is independent of W, we can see that the right-hand side of fIG.ll) 
is zero. Thus, ( 17. ip is obtained. 
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